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. We study the dynamics of r-modes in the ocean of a magnetic neutron star. 

On , We modeled the star's ocean with a spherical rotating thin shell and assumed 

. that the magnetic field symmetry axis is not aligned to the shell's spin axis. In 

1^ , the magnetohydrodynamic approximation, we calculate the frequency of £ = m r- 

O ■ modes in the shell of an incompressible fluid. Different r-modes with i and £±2 are 

^ . coupled by the inclined magnetic field. Kinematical secular effects for the motion 

O ! of a fluid element in the shell undergoing i = m = 2 r-mode are studied. The 

Q . magnetic corrected drift velocity of a given fluid element undergoing the i = m 

^ ! r-mode oscillations is obtained. The magnetic field increases the magnitude of 

^! the fluid drift produced by the r-mode oscillations. The drift velocity is strongly 

2 ! modulated by the inclined magnetic field. We show that the magnetic field is 

^ ', distorted by the high-£ magnetic r-modes more strongly than by the \ow-£ modes. 

! Further, due to the shear produced by the r-mode drift velocity, the high-£ modes 

• ^ ! in the ocean fluid will damp faster than the \ow-i ones. 

X' 
H : 

5t , Subject headings: stars: neutron - stars: magnetic - stars: rotation - stars: 

oscillations 



1. Introduction 

In a series of papers Bildsten et al. (1995, 1996) showed that neutron stars accreting at 
high accretion rates > 10~^ yr~^ are covered with massive oceans with densities > 10^ g 
cm~^. Both theoretical consideration of the nature of burning at these rates (Fushiki & Lamb 
1987; Bildsten 1993, 1995) and observational features of type I X-ray bursts (van der Klis 
1995) reveal that the star burns the accreting hydrogen and helium to mid-weight elements 
like C, O, Ne, etc., and accumulates a massive (~ 10~^ M©) degenerate liquid ocean. Fusion 
of this fuel to iron-group elements does not occur until densities > 10^ g cm~^. The neutron 
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stars with lower accreting rates (the atoll sources) burn the accreted hydrogen and helium 
directly to iron-group elements via the type I X-ray bursts which are crystallized by the 
larger Coulomb force at much lower densities. Further, they conjectured that the waves in 
these oceans might modulate the outgoing X-ray flux at frequencies comparable to what is 
observed. 

In this paper, motivated by Bildsten et al. conjecture we study the evolution of r- modes 
in the ocean of magnetic neutron stars. The r-modes those are analogous to the Rossby 
waves in the earth's oceans, are driven by the Coriolis force in the rotating stars. Their 
motions arc dominantly toroidal and their oscillation frequencies arc proportional to the 
rotation rate of the star, Q. Since even the small magnitude of Coriolis force drives r-modes 
in a rotating fluid, they should be considered in the dynamics of the fluid. We model the 
star's ocean with a thin rotating shell of incompressible inviscid fluid at the radius R which 
is sandwiched between two hard spheres. The role of spheres is to make sure that the fluid 
motion is restricted to a two-dimensional spherical surface. 

The paper is organized as follow. In the next section, we set up the magnetohydrody- 
namics equations for a spherical thin shell in the background of a uniform magnetic fleld 
in which magnetic fleld symmetry axis is not aligned with shell's spin axis. We flnd the 
magnetic coupling coefficients and then the corrected eigenfrequencies and eigenvectors for 
£ = m r-mode. In section 3 we consider the kinematical secular drift of r-modes, introduced 
by Rezzolla et al. (2000, 2001a,b), in a magnetic ocean fluid. We show that as well as the 
magnetic field strength causes the drift magnitude to increase, the inclination of the mag- 
netic field axis modulates the r-mode trajectories. The latter did not consider in the previous 
studies. Section 4 is devoted for further discussions. 



2. R-modes in magnetic rotating shell 

Historically, the study of oscillations of fiuids in a rotating shell in the presence of a 
magnetic field has been started about fifty years ago by geophysicists. Pioneering works 
by Hide (1966) and Stewartson (1966), are concerned with Earth's core magnetic modes 
that may be accounted to the observed secular changes in the main geomagnetic fleld at the 
Earth's surface. It is well known that the Earth's magnetic fleld at the surface drifts westward 
slowly, with the period of the order of one thousand years. They assumed an unrealistic 
magnetic field configuration characterized by a toroidal fleld of constant magnitude on the 
shell, and found the magnetic corrected toroidal motions might cause the field drift at the 
surface of the Earth. 
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In this paper we investigate the MHD perturbations of a uniformly rotating thin spher- 
ical shell of incompressible fluid of radius R, endowed with a uniform magnetic field (more 
realistic configuration) which is given by 

'B = Bp (cos 9 Sf — sin 9 e^) , (1) 

on the shell. Since the shell is sandwiched between two spherical hard covers, the motion 
of the fluid is restricted to a two-dimensional spherical surface. Further, we assume the 
symmetry axis of the fleld makes an angle /? to the shell's rotation axis. We work in an 
ideal MHD framework, so that the fleld lines are frozen into the fluid, and the magnetic fleld 
rotates at the same rate as the shell. 

For an uniformly rotating shell of incompressible fluid with angular frequency Vt around 
the z-axis, the fluid velocity perturbation fleld in the corotating frame is exactly determined 
by 

= a'i^^ , (2a) 
= ( -a^F,^(^,0)/sin^ e,- + a,F,^(^,0) ) e"--* , (2b) 

for any arbitrary function Uem{R) and with dispersion relation 

^^- = -777^- (2c) 



Here a is the amplitude of the mode. Equations (2) is the exact solution of the fluid 
equations of motion in fully nonlinear regime, see Levin & Ushomirsky (2001) for more 
detail and references therein. 

Recently, Morsink & Rezania (2002) (MR) have studied dynamics of a rotating star in 
the presence of a magnetic fleld in which its symmetric axis is not aligned to the star's spin 
axis. For the £ = m r-modes of an incompressible rotating fluid in the background of uniform 
magnetic fleld they obtained a shift proportional to the ratio of the magnetic energy to the 
rotational energy of the star in the r-mode eigenfrequencies, only. 

Following MR, we consider the magnetic force, produced by the fluid's perturbations 
of the rotating shell, drives the motion of fluid elements. As they've shown, the eigenvalue 
equation for eigenfrequencies and eigenvectors of the normal modes of fluid in the presence 
of magnetic field will be corrected by 

^a\ca - ^ ^ i^adCd \ = ctaCa, (3) 
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where the ca and cr^ are the magnetic corrected eigenvector and eigenfrequency in the 
rotating frame, respectively. Here is the energy of the mode in the rotating frame. The 
magnetic couphng coefficients, kad, can be calculated on the surface of sphere as 

i^AD = --r\-; [ SB*^ ■ 6Bd , (4) 



SBa = Vx (^^ X B) , (4a) 

where M — Bp/An is the magnetic field energy density ^. Note that in the above calculations 
we expand any perturbation vectors in terms of the modes of rotating star, as 

A 

Here upper-case Latin subscripts are used to label solutions and corresponds to the unique 
set of quantum numbers which describe the solutions. 

To calculate the magnetic couphng we use the method that is described by MR. They 
used the spin-weighted formalism in which: (a) the vectors are expressed in the orthonormal 
but complex basis {e+, e_} that is related to the usual orthonormal basis {eg, e^} by 

e± = (ee ± e^) ; (6) 

(b) all functions are expanded by spin-weighted spherical harmonics, ^£rnid,(t>), that are 
defined by (Campbell 1971) 



sYUe, 0) = ^ ^^^di.Me'""' (7) 

where the dl^{9) are the matrix representations for rotations through an angle 9 discussed 
in detail by Edmonds (1974). When s = 0, the spin- weighted spherical harmonics reduce to 
the regular spherical harmonics. The orthogonality relations are 



/ 



^i'm' yimd^ — Su'Smm'- (8) 

'An 

We use the convention that aYx^ii^, 0) = (-1)'"+'' -aYx-^,{9, 0). 

A uniform magnetic field tilted by an angle (3 from the shell's spin axis can be written 
in spin-weight formalism as 

m=l 

B{9, 0) = J] ( 6o'"(«) yifn{9, 0)eo + bf{a) ^,rn{9, (p)e- + h'fia) -Yirn{9, 0)e+ ) (9) 



m=— 1 



^Equation (1) leads V x B = 0. As a result, the term containing J vanishes. The term containing ^ x B 
does not vanish, but its contribution is smaller by ^I'^B? /(? for r-modes and we neglect it in equation 4. 
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where the functions bl''^{a) are given by 

b!f^^bl'^^^dUP), (10) 

and the d!^Q are given by 

c?io(/9) = ;^sin/?, cos/3, d\M^--^smp. (11) 

To rewrite equations (2) in spin-weighted formahsm, at first we normahze ^ such that the 
energy density of the modes in the corotating system is = T, where T is rotational kinetic 
energy density of the shelL Therefore the r-modes of non-magnetic incompressible thin shell 
have the form 

^niA^ no ^ 

""'""^^ ^ ~UiI+T) ^^^^^ 



'ti(G + i) 



ee^m^ =T= ^pR^Q' , (12c) 
The spin- weighted decomposition of the 74th r-mode is 

U^) = fi ^e,m,e_ + ^e,m,e+, (13) 
where the functions /^(r) are given by (Schenk et al. (2001)) 

Different components of the perturbed magnetic field, SB a, due to the r-modes in the thin 
shell, can be found in spin-weighted formalism as 



E^<%^a;> (15) 

C{iA,l,X)( ^ -)fn ^1 M V/^(l + (-l)^+^^), (15a) 

m— — 1 ^ ^ ^ ^ 

(15b) 



(-1) 

+1 



A,- 



X lA l\ i^^ r^^A ^ 

ruA m I a 



A + l) ( X (a 1\ {£A-mA + 2) / X Ia 1 
2 l-lOl/ V 2 1-12-1 



(15c) 
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( I k X 

where i ^ ) is a Wigner 3-j symbol (Edmonds 1974) and the constant 

C(£aJ,A) is defined by 



In equation (15) A takes on values Ia^'^a ± 1, although Eq. (15a) will survive only for 
X — Ia- Since the spin- weighted spherical harmonics obey the orthogonality relation (8), 
integrals over all angles of quantities quadratic in 5B have the simple form 



f^B^^^B^,* Sad + HB^'XiB^'X (l + (-l)'-**'")] , (ITa) 



where 5ad is the Kronecker delta function. The last step is due to the symmetry property 
(15b). It then follows that the magnetic coupling coefficients between r-modes is zero unless 
both modes have the same parity. Since the triangle inequalities — 1<A<£^-|-1 and 
^D — 1<A<^£)-|-1 must be satisfied, only modes satisfying Id — (-aAa ± 2 have non- 
zero magnetic coupling. Evaluating the terms appearing in equations (15a) and (15c) for 
I a — in A, the individual terms for the allowed values of A are 

m=— 1 \ \ / X 



V TO= — 1 ^ 



{I A + fh) I A fh 



(18d) 



The magnetic coupling coefficients between different r-modes can now be computed. A 
straight-forward calculation yields the self-coupling term, kaa as 

>^AA = - 4(2^^"+^) {^^(^^ + i)(£a + 2) + 3 + ^ [^A{^A + 1)(2£^ + 2^A - a) - 9] sin^ /?! . 

(19) 



- 7- 



The off-diagonal term, k,a,a+2 and ha,a-2 will be 



...... = = + + 3)^ sin^ . (20a) 



Therefore, the equation of motion for mode's expansion coefficients of the magnetically 
modified >lth r-modes reduces to 



(^aCa — ^A 



Ka,A-2\CA-2 + 1 + -^Z^ll^AAl CA + ^\i^A,A+2\CA+2 



(21) 



The final equations (19)-(21) give the magnetically corrected eigenfrequencies a a and eigen- 
vectors ca for the r-modes. To compute these eigenfrequencies and eigenvectors one should 
solve the N dimensional linear algebraic system (21), where N is the number of the modes. 
Therefore, for a given = 1, 2, • • •, the N x N matrix should be diagonalized. Since it is not 
possible to solve an infinite dimensions system in practice, we truncate the infinite matrix 
to a finite matrix and large enough N. Further we note that for the lower inclination angles, 
the magnetic coupling coefficients \kad \ increase by £% while for the larger inclination angles 
by £\. Since we are working in the linear regime, we have to keep {M./2T)\k,aa\ < 1- The 
ratio of magnetic energy density to rotational energy density in the covered ocean at the 
surface of star is 

where p is the ocean's density. In tables 1 and 2 we present the eigenfrequencies of £a — 
niA — 1, 2, • • • , 6 r-modes. In table 1 we assume the inclination angle between field axis 
and shell's spin axis is zero, P — 0, and calculate the eigenfrequencies cr^ for different ratio 
of the magnetic energy density to the rotational energy density, A1/2T, by increasing the 
magnetic field strength, see equation (22). In table 2 we compute the eigenfrequencies by 
increasing the inclination angle for a specific energy ratio M./2T = 10^^. As a result, the 
eigenfrequency of the mode increases by increasing both magnetic field energy density and 
inclination angle. 



In the next section we consider the kinematical secular effects of r-modes in the oceanic 
fluid. 
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3. The r-mode drift by the magnetic field 

In a series of papers Rezzolla et al. (2000, 2001a,b) showed the existence of kinematical 
secular velocity field of r-mode oscillations which interacts with the background magnetic 
field of the star. Although, this interaction is non-linear and should be considered in non- 
linear regime, one can find some second-order quantities from linear result (Rezzolla et al. 
2000). As an example, the second-order secular toroidal drift which appear on isobaric 
surfaces, can be obtained through the linear velocity field. 

To investigate this secular effect in the fiuid motion undergoing the magnetic corrected 
r-mode oscillations in the ocean fiuid, we consider the real part of the magnetic corrected 
fiuid velocity of — in a r-mode, equations (2), 

e{t) = - ^ a g(£^) ca a a (sin OY^'^ cos(£a0 - 271(7^^), (23a) 

A 

4){t) = ^ Oi Q{£a) ca aA cose {smef^"^ sm{eA(p - 27raAt), (23b) 

A 

where a a = cfa/^ and 

Q{^a) = (^-0 ■^V^£a(^a + 1)(2£a + 1)!/4. (23c) 

In equations (23) a a and ca are the magnetic corrected eigenfrequency and eigenvector, 
respectively, and the time, t, is in unit of the star's period, P. 

In figures 1 to 3 we plotted the result of numerical integration of equations (23) in the 
shell's northern hemisphere for the magnetic corrected I a — n^A — 2 r-mode with amplitude 
a = 0.1. Note, we assumed that at t = the only £a = tua — 2 r-mode is present in the 
fiuid and all calculations are done in the corotating frame. 

In the case of zero inclination angle, /3 = 0, our results for all mode's latitudes are in 
good agreement to those results which have been reported by Rezzolla et al. (2000, 2001a). 
In figure 1 we plotted 9{t) vs (j){t) from t = to t = 7.45P (almost 5 cycles of the mode), 
with initial values 0(0) = 0, ^^(0) = 7r/2, I a = f^A = 2, zero inclination angle /5 = 0, and 
two different values oi M./2T = and 10~^ The aligned field causes the drift in 0-direction 
(produced by the r-mode) to increase. The corrected expression for the drift velocity of a 
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given fluid element in the background of the uniform magnetic field will be ^ 



(24) 



A 



KA{d) = Q\tA)c\{sme) 



sm^e -2{Ia-1)cos^0 



(24a) 



As a result, because of nonzero inclination angle, the drift velocity of a fluid element under- 
going the Ia r-mode oscillation is affected by all ± 2 r-mode oscillations. Equation (24) 
reduces to one obtained by RezzoUa et al. (2000, 2001a), for the non-magnetic I a = "m-A = 2 
r-mode, ie. = K2{0)a(t)uj2{t)etf,. Therefore, the total displacement in 4> from the onset of 
the oscillation at tg to time t is 



The change in due to the magnetic field is given in tables 3 and 4. As a result, A(f) = 



magnetic r-modes are drifted more than the \ow-£a ones. 

As shown in figures 2 and 3, the r-mode trajectories in both 9- and 0-directions are 
strongly modulated by the inclined magnetic field. In figure 2 from left to right, we plotted 
time evolution of r-mode displacements in ^-direction, in ^-direction, and the projected 
{9, (f)) r-mode trajectory, respectively. All panels correspond to A4/2T — 10~^ and are 
plotted with initial conditions (f)(0) — 0, ^(0) — 7r/2. The inchnation angle is P — for 
the top, P = n/6 for the middle, and P = n/2 for the bottom panels. The larger the 
field's angle the stronger the modulation. First and second columns of figure 2 show how 
9- and (p- displacements of the initial Ia = ^t^a = 2 r-mode are affected by the inclined 
field, respectively. In the third column we show the motion of fluid element in the northern 
hemisphere of the shell undergoing the £a = f^A = 2 r-mode oscillations. The projected 
trajectories 9{t) sin 9{t) cos (f){t) and (pit) sin 9{t) cos 0(t) show how the inchned magnetic field 
change the motions which will be detected by a corotating observer in the rotation equator 
of the shell. 

It is interesting to note that the results show the modulations are more sensitive to the 
field inclination angle rather than the magnetic field strength. This can be understood from 
equations (23). A non-diagonal matrix of the corrected eigenvectors ca causes modulations 
in both 9 and (f) trajectories. For any non-zero value olM./2T ^ 0, no matter how small, and 




(25) 



mag 



nonmag 



increases by increasing the magnetic energy density. Further the high-^^ 



^Following RezzoUa et al. (2000, 2001a) equation (24) is obtained by expanding equations (23) in powers 
of a, averaging over a gyration, and retaining only the lowest-order nonvanishing term. 



-10- 



/3 7^ the eigenvectors matrix is not diagonal. Furthermore, the magnitude of the corrected 
eigenvector ca changes significantly by varying while it is almost constant for different 
(and non-zero) values of hAjTI . This fact is shown in figure 4. In the left panel, we plotted 
the element C2(2, 2) as a function of /? for a fixed value of M./2T ^ 10~^*^. As is shown, by 
increasing j3 the value of C2(2,2) drops by ~ .02. While, the change in the element C2(2,2) 
due to the M./2T for a fixed value /5 = 7r/2 is < 10~®, which is shown in the right panel 
of figure 4. Therefore, no matter how the value of the energy ratio small, the fluid motions 
are strongly modulated by a magnetic field making a large inclination angle to the star's 
symmetry axis. 

Furthermore, the larger the fiuid element's latitude, the stronger the modulation. This 
feature is shown in figure 3 which is plotted by different initial conditions 0(0) =0, 6'(0) = 
7r/6 for both /3 = (the top panels) and (3 — 7^/2 (the bottom panels). Comparing both 
bottom panels in figures 2 and 3 reveals that the magnetic field (having the same energy 
density and inclination angle) has less effect on the lower latitude modes. 

Finally we note that A0 also increases by increasing the inchnation angle, see the third 
column of figures 2 and 3. Therefore, the motion of a fluid element of the shell undergoing 
the magnetic r-mode oscillation is drifted more than nonmagnetic r-mode oscillation both 

by field energy density and field inclination angle. These properties can be understood by 
considering both and are proportional to the corrected eigenfrequency and eigenvector, 
see equations (24) and (25). The corrected eigenvector changes by changing inclination angle 
while the corrected eigenfrequency increases by increasing field energy density. 

The toroidal fiuid motions produced by r-mode are perpendicular to the magnetic field. 
They distort the field and increase the field energy density. As shown by Rezzolla et al. 
(2000) if the energy that it needs to distort the magnetic field required by the mode is 
greater than the mode energy during an r-mode oscillation, the distortion of the magnetic 
field which is proportional to the v^, will prevent the oscillation from occurring. To make an 
estimate at which field strength the distortion effect might get to be important we calculate 
the ratio of the change in the magnetic energy density of the shell, 5E^ — J (55^/47r)dQ, to 
the mode energy density during one oscillation as 

SE^ M 



, 

eA 2T 



\haa\ , (26) 



^^f B V /lO^ gcm-3\ /lO km\ VlO^ HzV r £3 for small /3 , 
" \ -p )\-^) \-^) in for large/? ^2'^) 

Thus 5E^ > if S > Senticai ~ 3 x p^^ Rw^^(\, where 77 = -3/2 and 77 = -2 for 
small and large inclination angles, respectively. Here pg = p/lO^ g cm~^, Riq — R/10 km, 
and Qs — Q/IO^ s~^. Hence the smaller field strength would prevent the high-f^ r-mode 
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oscillations. Further, equation (26) shows that for the higher field inchnation angle the 
distortion effect will be stronger. 

In addition the increase in the energy of the magnetic field produced by the r-mode drift 
reduces the mode energy and causes damping. The rate that the magnetic field extracts 
energy density from the mode can be estimated by 

(l/47r) / B^{dB^/dt)dn « KA{e)a^aABpB^ . (27) 

The azimuthal field B^ is generated by r-mode drift velocity from the background poloidal 
field. The change in can be calculated approximately for the magnetic r-modcs from 
the time that the oscillation begins to time t by AB^p^t) ^ J^^^ KA{0)a'^{t)aA{t)Bpdt 
(RezzoUa et al. 2000). Equation (27) shows that the energy loss from the mode increases 
when either the field energy density or the inclination angle is increased. Furthermore, since 
the high-^^ magnetic r-modes are drifted more than the \ow-£a ones, they will damp faster 
than the low-£^ modes. 



4. Discussion 

In this paper, by assuming neutron stars in the low-mass X-ray binaries (LMXBs) 
are covered by degenerate massive oceans, we have studied the interaction of the inclined 
magnetic field with r-mode oscillations of an incompressible ocean fiuid in a thin shell ap- 
proximation. Our analysis shows that in the case of large inclination angle, even for weak 
magnetic field energy density in comparison with the star's rotational energy density, £ and 
£±2 r-modes are coupled significantly. 

The kinematical secular drift produced by the r-mode, is studied in the background of 
the inclined magnetic field. The results show that the magnetic field tends to increase the 
fiuid drift in the 0-dircction. As a result, for the same time period, the fluids in the presence 
of the magnetic field undergo larger drift than fiuids with no background magnetic field. 

Since the magnetic field is distorted by the r-mode, we showed that the field distortion 
by the high-£ magnetic r-mode is stronger than the low-^ modes. Thus a smaller critical 
field is needed to prevent the high-£ r-mode oscillation from occurring. Further, the shear 
produced by the drift velocity of the magnetic r-mode will damp high-£ modes faster than 
the \ow-i ones. Note that in the present paper, we assume the r-mode is stable, so we do 
not deal with its unstable growth in the magnetic star as has been considered by RezzoUa 
et al. (2000). 



- 12 - 



Observational and phenomenological studies of more than twenty LMXBs reveal that 
these accreting neutron stars arc source of millisecond oscillations, burst oscillations, and 
quasi-periodic oscillations in Hz-kHz range. Different authors proposed different possible 
mechanisms to explain the observations in LMXBs, but the regularity of these oscillations 
in many sources led the investigators to consider radial or nonradial oscillations in a coronal 
accretion flow near the Eddington limit, see the review by van der Klis (2000). By assuming 
neutron stars in LMXBs covered by an ocean fluid Bildsten et al. (1995, 1996) showed 
that both low and high radial orders of ocean g-modes might be responsible for observed low 
frequency (~ Hz) oscillations in these stars. Recently Heyl (2001) discussed different possible 
oscillations (the g-modes, the Kelvin modes, and the buoyant r-modes) that may excite in 
the neutron star ocean during a Type-I burst, and concluded that the buoyant r-modes may 
meet the observations of low frequency oscillations in Type-I X-ray bursts quite well. 

Figures 2 and 3 show how the inclined magnetic fleld modulates iA — tua — 2 r-modes. 

The associated envelope frequency, /envelope, decreases as /envelope ~z^s/10 — z/^/lG when the 
inclination angle increases from ~ 7r/6 — 7r/2. Here Vs = ^/2it is the star's spin frequency 
. The frequencies of these new modes are comparable to those observed oscillations in the 
LMXBs in low frequencies. For the star's spin frequency Ug ~ 300 Hz, the wave envelope 
frequency is /envelope ^ 20 — 30 Hz. Different frequencies may be obtained by different 
initial Ia = tua r-modes. Further, the amplitude of the modulated oscillations increases by 
increasing the inclination angle and/or mode's altitude. 

Finally wc note that during the evolution of the r-mode instability, one would expect the 
growth of the amplitude of the new modes. This feature can be understood from equations 
(23) which show that both 6 and trajectories of the r-mode are proportional to the mode's 
amplitude, a. Further, changes in the star's spin frequency in this period, lead to change 
in the frequency of modulated modes as well. In flgure 5 we plotted time evolution of 0{t)- 
displacement of the Ia — tua = 2 r-mode, from t = Q to t — lOOP, for different initial 
mode's amplitudes a = 0.1, 0.5 (the top panels) and a = 0.8, 0.9 (the bottom panels). 
All panels corresponds to M./2T = 10~^, f3 = 7r/2 and are plotted with initial conditions 
0(0) = and ^^(0) = 7r/2. By increasing a both amplitude and period of the modulated 
oscillations increase. Further, for a close to the saturation limit, ttgat = 1-0, the behavior of 
the modulated oscillations changes drastically. Therefore, to analyze the problem carefully, 
one needs to solve both the nonlinear equations (23) and the equations governing the r-mode 
evolution together. This is currently under investigation (Rezania 2002). 

Though there is still no firm observational evidence for the low magnetic field strength 
in the LMXBs, these objects are beheved to be neutron stars with weak magnetic fields 
{B < 10^ G), high accretion rates (M > 10~^ M0 yr~^), and milliseconds periods (van 
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der Klis 2000). Our calculations show the existence of an inclined (even weak) magnetic 
field modulates the normal modes of an ocean fluid to a new mode with frequency /envelope- 
Therefore, by observing normal modes in the neutron star's ocean, one would expect to 
observe these new modes. 

The author is grateful to S. M. Morsink, S. Sengupta and L. RezzeoUa for reading the 
manuscript and numerous discussions. I would like to thank referee for his useful comments. 
This research was supported by the Natural Sciences and Engineering Research Council of 
Canada. 
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Table 1. The eigenfrequency of ^A = f^A r-modes in the thin shell in the presence of a 
uniform magnetic field, a^/fi. The angle between magnetic field axis and shell's spin axis 
is (3 = 0. The energy in the magnetic field is A4 = Bp/An and the rotational energy is 
T = pa'^Q? 12. A number a x 10=*=^ is written as o ± 6. 



M/2r 


^A = l 


^A = 2 


^A = 3 


£^ = 4 


^^ = 5 


^^ = 6 


0.0 


-.10000+01 


-.66667+00 


-.50000+00 


-.40000+00 


-.33333+00 


-.28571+00 


1.0-04 


-.10000+01 


-.66686+00 


-.50004+00 


-.40056+00 


-.33415+00 


-.28684+00 


1.0-03 


-.10005+01 


-.66943+00 


-.50838+00 


-.41797+00 


-.37886+00 




1.0-02 


-.10090+01 


-.68595+00 


-.53500+00 


-.45591+00 
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Table 2. The eigenfrequency of ^A = r-modes in the thin shell in the presence of a 
uniform magnetic field, (7^/5^ for different value of inclination angle, /3. The ratio of 
magnetic energy to the rotational energy is A4/2T — 10~^. A number a x 10='='' is written 

as a ± 6. 





£a^1 


£a^2 




0.0 


-.10090+01 


-.68595+00 


-.53500+00 


7r/6 


-.10080+01 


-.68874+00 


-.54775+00 


7r/4 


-.10068+01 


-.69038+00 


-.55827+00 


7r/3 


-.10056+01 


-.69155+00 


-.56837+00 


7r/2 


-.10043+01 


-.69250+00 


-.57834+00 



Tabic 3. Change in displacement 0, 

= (pmag " (t>nonmag for £a = "^A = 2 r-Hlode, due tO 

the magnetic field, from t — Q to t — TAhP. Initial 
values 0(0) = 0, ^(0) = 7r/2 and zero inclination angle 
/3 = are considered. A number a x 10='='' is written as 

a±h. 



M/2T A(f) 



0.0 0.0 

1.0-04 .36403-06 

1.0-03 .19299-05 

1.0-02 .29707-03 

1.0-01 .10390-01 
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Table 4. Change in displacement 0, 
= (pmag " (pnonrnag for different ^A = tua r-mode, 
due to the magnetic field, from t = to t = 7.45P. 
Initial values 0(0) = 0, ^(0) = 7r/2, M/2T = 10'^ and 
zero inclination angle /3 = are considered. A number 
a X 10='='' is written as a ± 6. 



1 .56340-06 

2 .29707-04 

3 .39012-03 

4 .14556-02 
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Fig. 1. — Motion of a fiducial fluid element in the northern hemisphere of a shell undergoing 
iA — TTT'A — 2 r-mode with amphtude a — 0.1. The projected trajectories are plotted for 7.45 
of period of star (almost 5 cycles of the mode), with initial values 0(0) = 0, 9{0) = 7i/2 and 
zero inclination angle (3 = 0. The dash sign refers to zero magnetic field, Ai = 0, while the 
continuous line corresponds to assumed value A4/2T = 10~^. Both 9 and are in radian. 
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Fig. 2. — Motion of a fiducial fluid element in the northern hemisphere of a shell undergoing 
= "m-A = 2 r-mode with amplitude a = 0.1. The projected trajectories are plotted for 20 
of period of star (almost 13 cycles of the mode), with initial values 0(0) = 0, ^(0) — 7r/2 
and M./2T — 10"^. Both 9 and are in radian and t is in unit of star's period, P. (a) The 
top panels are plotted for P = 0, (b) The middle panels are plotted for /3 = 7r/6, (c) The 
bottom panels are plotted for = 7r/2. 
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Fig. 3. — Same as figure 2 with initial values 0(0) = 0, ^(0) = 7r/6. (a) The top panels are 
plotted for /3 = 0, (b) The bottom panels are plotted for /3 = 7r/2. 
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C2(2,2) 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



/3 



C2(2,2) 



+ + 



1e-22 1e-18 1e-14 1e-10 1e-06 1e-2 

M/2T 



Fig. 4. — Change in the element 02(2, 2) as a function of the inchnation angle (5 and/or the 
energy ratio M/2r. (a)The left panel is plotted for M/2r ^ 10'^°. The values of /3 in 
radian, (b) The right panel is plotted for P — 7r/2. Both axes have logarithmic scale. 
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Fig. 5. — Time evolution of 6{t) — 6{0) displacement of the = = 2 r-mode, from t = 
to t — lOOP, for different initial mode's amplitudes a. (a) The left-top panel is plotted for 
a — 0.1, (b) The right-top panel is plotted for a — 0.5, (c) The left-bottom panel is plotted 
for a — 0.8, (d) The right-bottom panel is plotted for a — 0.9. All panels correspond to 
M/2r = 10-^ (5 = 7r/2 and are plotted with initial conditions 0(0) = and ^(0) = 7r/2. 



